Non-linear metric perturbation enhancement of primordial gravitational waves 



O 

-)— » 
o 

O 



M. Bastero-Gil/^2,g j^ Macias-Perez,2'[t| and D. Santos^'S 

' Departamento de Fisica Teorica y del Cosmos, Universidad de Granada, Granada- 18071, Spain 

^Lahoratoire de Physique Subatomique et Cosmologie, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, 

Institut Polytechnique de Grenoble, Grenoble-38026, France 

We present the evolution of the full set of Einstein equations during preheating after inflation. We 
study a generic supersymmetric model of hybrid inflation, integrating fields and metric fluctuations 
in a 3-diniensional lattice. We take initial conditions consistent with Eintein's constraint equations. 
The induced preheating of the metric fluctuations is not large enough to backreact onto the fields, but 
preheating of the scalar modes does affect the evolution of vector and tensor modes, fn particular, 
they do enhance the induced stochastic background of gravitational waves during preheating, giving 
an energy density in general an order of magnitude larger than that obtained by evolving the tensors 
fluctuations in an homogeneous background metric. This enhancement can improve the expectations 
for detection by planned gravitational waves observatories. 
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Cosmic microwave background (CMB) measurements 
[l| , are consistent with an early period of inflation, which 
gives rise to the primordial curvature perturbation which 
seeds the large scale structure observed today. The ob- 
served power spectrum of temperature anisotropics is 
consistent with a Gaussian, adiabatic primordial spec- 
trum, and present data sets at most an upper limit on the 
level of the primordial tensor contribution predicted by 
inflation. This contribution could be observed through 
B-mode polarization by the European Space Agency's 
Planck mission [2] and future CMB polarization experi- 
ments. Inflation should be followed by a reheating period, 
during which the inflationary vacuum energy is converted 
into radiation. During the first stages of reheating, the 
evolution of the system may be dominated by nonper- 
turbative effects such as those of preheating, i.e., para- 
metric amplification of quantum field fiuctuations in a 
background of oscillating fields [3|. Through paramet- 
ric resonance, field mode amplitudes grow exponentially 
within certain resonance bands in k space, this being an 
efficient way of transferring vacuum energy into radia- 
tion ;3]. It does also enhance the tensor perturbations, 
sourced by the field anisotropic stress-energy tensor, giv- 
ing rise to a stochastic background of gravitational waves 
(GWs) [i^. 

In most of the studies of preheating of GW, fields 
and tensor fluctuations are evolved in a background 
Friedmann-Robertson-Walker (FRW) metric. However, 
beyond linear perturbation, tensors are seeded by the 
other components of the metric, scalar and vectors Q- 
Given that the metric and fields are nonlinearly coupled 
through the Einstein equations, the parametric ampli- 
fication of field fiuctuations will be rapidly transferred 



to all metric perturbations. In a previous paper j8| we 
showed that this is indeed the case for the scalar metric 
perturbations. We showed that during the resonance the 
scalar metric source term for the tensors could be compa- 
rable to the anisotropic field stress tensor, arguing that 
they could affect the amplification of tensors. The pur- 
pose of this Letter is to show that this is indeed the case 
for hybrid models independently of model parameters. 
For that aim we have integrated the full set of Einstein 
equations, with the stress-energy tensor provided by the 
fields in hybrid inflation. 

Einstein equations are written in the so-called 
Baumgarte-Shapiro-Shibata-Nakamura formalism [9|, 
widely used because it provides a stable numerical evolu- 
tion. One starts with the Arnowitt-Deser-Misner metric, 
given in terms of the lapse TV and shift vector N^ func- 
tions, and the spatial metric jij-. 



ds^ = N^dt^ - -i^J{dx' + N'dt){dx^ + N^dt) , 
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and performs a conformal transformation of the spatial 
metric, ^ij = exp{—4:l3)"fij , with detjij — 1. The metric 
dynamical variables are then /? and jij , the trace of the 
extrinsic curvature K — j^^Kij, and its traceless part 
Aij = exp{—4l3){Kij'-jijK/i). In addition, a connection 
variable I"* — T^f^j-''^ = — 9j7-'* is introduced in order 
to compute the Ricci curvature more accurately. The 
lapse and shift vector are gauge functions, and we choose 
to work in the synchronous gauge and set A^ = 1 and 
N' = 0. In this gauge we have K = 6/3, with {K)/3 = 
H{t) being the average expansion rate, and {e^^) = a{t) 
the average scale factor ( "(•)" denotes spatial average.). 
The spatial metric jij encodes the two transverse and 
traceless degrees of freedom of the gravity waves, plus one 
additional scalar mode and 2 vector degrees of freedom, 
with 7jj- = 2 Aij. 

The matter source is given by the stress-energy ten- 
sor of the two scalar fields (which we take as real scalar 
fields) in supersymmetric hybrid inflation, with potential 



^(0, X) = Fo +g2;^V4 + g2($2 _ ^2)^2 ^ ^2^2/2, where 
$ is the inflaton field and x the waterfall field which trig- 
gers the phase transition at the end of inflation once the 
inflaton goes below the critical value </>£. The vacuum 
energy Vo is adjusted such that the potential vanish at 
the global minimum, Vq — g'^cf)^, and the mass for the in- 
flaton field is set by the choice of the slow-roll parameter 
77^ = rn^m/^/Vo, where mp is the reduced Planck mass. 
The set of equations of motion are then given by the 
Klein-Gordon equations for the fields in the curved space 
given by the metric ([T]) , the equations of motion for the 
connection F* , and those for the metric variables Q given 
by 
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where T — {9^ -I- x^)/2 is the kinetic energy of the 
fields; the field-dependent source is Mij = mp^{di^dj^+ 
dixdjx), Rij is the Ricci tensor of the metric 7.^ , and the 
superscript "TF" denotes the trace-free part of the ten- 
sor. 

The system is placed in a finite and discrete 3D comov- 
ing box of length L and A'^ sites per spatial dimension. 
The procedure introduces a comoving ultraviolet cut-off 
in both space and momentum. Ideally, one would like to 
keep both of them large enough: A big comoving box, 
larger than the initial Hubble size L > H{0)^^, ensures 
that the observable initial Universe is within the box; 
but a large enough number of points is also required to 
have all the relevant momentum modes up to 0{N/L). 
For the problem of preheating after inflation, one tries 
to optimize the choice of this ratio 0{N/L) to have an 
ultraviolet comoving momentum cut-off still larger than 
the preheating cut-off by the end of the simulation. This 
means that already at the start of the simulation our 
comoving box is smaller than the observable Universe, 
but the relevant physical modes for preheating are all 
included. 

For the fields we follow the standard procedure, with 
the quantum field theory being replaced by an equiva- 
lent classical field theory. The quantum origin of the 
fields remains in the stochastic nature of the initial con- 
ditions^. We start the simulations some fraction of e- 
fold AA'e — 0.05 after the end of infiation, with the 
background inflaton field still close to the critical point, 
($) = (/)cexp(— ry^AiVe), and its background velocity 
given by the slow- roll conditions^. The background val- 



A description of how to follow the quantuin-to-classical transition 

for the scalar fields, and how to fix the initial conditions for the 

modes in the classical regime can be found in [10]. Nevertheless, 

we have checked that for the particular problem at hand, a better 

determination of the initial conditions has little impact on the 

results. 

For the model parameters considered, the inflaton field reaches 
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FIG. 1: Evolution of the energy density of GWs, paw, nor- 
malized by the initial vacuum energy pi, for different values 
of the box size. Solid lines are the result of integrating the 
Einstein equations; dashed lines show the evolution of paw 
without including other metric fluctuations. 



ues for the waterfall field are set to zero at this point. 
Classical inflaton field fluctuation in a spatial lattice with 
periodic boundary conditions are expanded as usual in 
Fourier modes $/£, with an initial vacuum amplitude 

\M0)\ 



l/^/2uJk, where wj, = y/k^ + mA. The re- 
maining variables must be chosen to satisfy the Einstein 
constraint equations, the momentum and the Hamilto- 
nian constraint, at t = 0. This is the well known initial- 
value problem in general relativity [ll|. We choose ini- 
tially vanishing tensors and vectors, i.e., 7ij(0) = Sij, 
such that the momentum constraint reduces to 2diK — 
— {i / m?p){f^di^ -\- xdix) ■ In order to fulfill this equation, 
we set the waterfaU field as x(0) = $(0) - (<&(0)), and^ 
x(0) = ($(0)) - <i)(0). This allows us to solve the mo- 
mentum constraint for the initial value of the expansion 
rate fluctuations, and use the hamiltonian constraint to 
flx the fluctuations of the scale factor^. 

Preheating in hybrid inflationjmodels has been exten- 
sively studied in the literature [131] . The parametric am- 
plification of the fluctuations takes place first through an 



the critical point still in the slow-regime, and in the absence 
of any extra mechanism that could speed up the inflaton, slow- 
roll is a good approximation across the bifurcation point. For 
the dependency of the GW spectra on the initial velocity of the 
inflation fleld see Qj. 

Given that once preheating starts the inflaton and waterfield fluc- 
tuations become comparable, results are independent of which 
field we set initially in the vacuum, and which is set by the con- 
straint. 
* Violations of the constraints arc unavoidable during numerical 
integration, and their growth makes the system unstable. Nu- 
merically, we have checked that adding the covariant derivative 
of the momentum constraint in Eq. ((21 renders the system more 
stable ^IJ]. 



spinodal instability for the fields, during which the lower 
modes are quickly amplified. After a few oscillations, the 
amplitude of the fields has decayed enough to be out of 
the spinodal region, and tachyonic preheating ends. We 
just follow the evolution of the fields and metric vari- 
ables up to the end of the resonance, before we lose the 
ultraviolet cut-off for the field modes. Metric variables 
will follow the same pattern of parametric amplification 



as the fields, with 



rup |$| . Keeping only the 



leading /? terms in the Ricci tensor, we have: 
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and thus RJa^ becomes comparable to the field contribu- 
tion Mjj^ in Eq. ©. 

The metric variable Aij is traceless but nontransverse, FIG. 2: Comparison of the evolution of the absolute values 

i.e., it contains more degrees of freedom than those two of the source terms Mj/ and R^^y in Eq. (g]), in logarithmic 

corresponding to GWs. The traceless and transverse scale. The spikes in the plot are localized where the corre- 

(TT) components are projected by using the operator sponding oscillating term vanishes. 
i Ay,,„(fc) = Pa{k)Pjm{k) - P,,{k)Pim{k)/2, where 
Pij = Sij — kikj and ki = ki/k. The energy density 



of the GW is then paw = m%{A,j{t,x)A'i{t,x))^^ = 
mlJd^k\Af/'{t,k)\^ 0, i]. In Fig. [I] we have plotted 
PGWj normalized to the initial vacuum energy pi = g'^'f't- 
We have taken as parameter models g — 0.01, (pc — 
0.005mp and ry^ = 0.05. The value of pcw does not de- 
pend on the value of the coupling, which can be rescaled 
out from all the equations, but it does depend on the 
value of cj)c [J-[a| , which sets the scale for the field source 
term^ such that paw/ Pi oc <j)1- We have included the 
results for different choices of the box size, to show that 
the final value of pcw does not depend on the choice of 
L as far as we have all the relevant modes to start and 
end tachyonic preheating. With a larger comoving box 
we have more modes in the low momentum regime, and 
then tachyonic resonance starts slightly sooner, as can be 
seen in the plot. 

We compare in Fig. [T] the results obtained when inte- 
grating the full Einstein equations (solid lines) with those 
obtained when integrating the tensor modes in a FRW 
background metric (dashed lines) . The former are always 
roughly an order of magnitude larger due mainly to the 
contribution of the scalar modes of the metric fluctua- 
tions in Eq. (|4]). In Fig. [2l we show that both source 
terms, M^^ and the leading term Rf^^ in Eq. (0]), are of 
the same order by the end of the resonance. Scalar met- 
ric fluctuations start growing immediately after inflation 
due to the increase in the kinetic energy of the held [see 
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solid: with P{x), yM _ 
dashed: with a(t) 




The parameter 0c also sets the initial value of the expansion 
rate Hi/{g<fic) = <t>c/{V3mp), and relatively large values of Hi 
have an impact on the evolution of the resonance, due to the fast 
redshifting of the field amplitudes. Tachyonic resonance ends 
once these amplitudes go below a certain value, and therefore 
the larger (f>c and Hi , the shorter resonance time, and the smaller 
the amplification of the field fluctuations and, thus, of GWs. 
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FIG. 3: Spectrum pawik)/ pi per logarithmic k interval, for 
{g(j>c)t = 60 to 200 in intervals of 20 from bottom to top. 
Solid lines are the result of integrating the Einstein equations; 
dashed lines show the result of the integration in a FRW met- 
ric. 



Eq. ([2])], and this effect leads the initial growth of the 
tensor perturbations. 

In Fig. [31 we show the spectrum of GW per logarith- 
mic frequency interval, PGw{k) = dpaw/dlnk, at dif- 
ferent times until the end of the tachyonic resonance. 
Compared with the spectrum obtained without metric 
perturbations, the enhancement in the spectrum is local- 
ized at around the peak of the spectrum, shifted towards 
lower values at around k ~ g(j)c. This effect at low mo- 
menta is again due to the scalar metric perturbations, 
which spectrum peaks below gcjjc- This can be seen in 
Fig. m where we compare the power spectra of the in- 
flaton field fluctuations P^ = |<I>(A:)p, with those of the 
metric fluctuations /3 and Af^ , at the end of the res- 
onance. The held spectrum at i = is included as a 
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FIG. 4: Comparison of the power spectrum for the infla- 
ton field P^{k), and metric perturbations Pis{k) and Pa, at 
{g4>c)t = 100 (lower curves) and {g4>c)t = 200 (upper curves). 
We have rescaled Pa by (f)'^''^ to have all the spectra in the 
same scale. For P/j, we have k^'^Pp 
turn constraint at t = 81. 



reference, to check that by the end of the simulation the 
highest momentum modes have not yet been amplified. 

In summary, we have shown that nonlinear effects due 
to metric perturbations enhance the amplitude of the 
GW stochastic background by an order of magnitude 
with respect to the calculations in a FRW background 
metric. Taking into account that PGw(k) is redshifted 
like radiation after the resonance, and assuming that en- 
tropy is conserved from reheating onwards, its maximum 
present-day value normalized by the critical density to- 
day is given by: 



h'nP^^ ~ 5.5 X 10- 



(t>c/mp 
0.005 



Pi 
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where Trh is the reheating temperature at which the 
Universe becomes radiation-dominated. This ampli- 
tude is within the reach of the future GW observa- 
tory Advanced-LIGO for (p^ > O.OObmp, or BBO for 
0c ^ 0.005r7ip jl4|. However, today's values for the 

frequency are / = 6.4 x lO^''^[k/{gq},)]{TRH/p]^^)^^^ 
Hz, while the operating frequency range for Advanced 
LIGO is 1 — 10'^ Hz, and BBO will operate in the range 
lO^'^ — 10^ Hz. Thus, to have the GW spectrum within 
the observable range, the coupling should be at most of 
the order of g ~ 10~^^. The numerical simulations so far 
can resolve the spectrum for the typical frequency of the 
resonance, k ~ gcjjc and above, but not the infrared tail 
for subhorizon and supcrhorizon modes. The behavior of 
the GW spectrum in this infrared range and how they 
are affected by nonlinear metric effects is still an open 
question. A tail of subhorizon modes at the time of pre- 
heating rising slower than k^ could be detected by BBO, 
although fpeak wcrc in the range of 10^ Hz. Similarly, in 
models where t he p rimordial curvature is enhanced dur- 
ing preheating Il5!|, one could expect a similar effect for 
the tensor primordial spectrum, up to the level detectable 
in CMB polarization experiments. 
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